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The fermionic Hubbard model plays a fundamental role in the description of 
strongly correlated materials. Here we report on the realization of this Hamil- 
tonian using a repulsively interacting spin mixture of ultracold 40 K atoms in 
a 3D optical lattice. We have implemented a new method to directly measure 
the compressibility of the quantum gas in the trap using in-situ imaging and 
independent control of external confinement and lattice depth. Together with 
a comparison to ab-initio Dynamical Mean Field Theory calculations, we show 
how the system evolves for increasing confinement from a compressible dilute 
metal over a strongly-interacting Fermi liquid into a band insulating state. 
For strong interactions, we find evidence for an emergent incompressible Mott 
insulating phase. 
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Introduction. Interacting fermions in periodic potentials lie at the heart of modern condensed 
matter physics, presenting some of the most challenging problems to quantum many-body the- 
ory. A prominent example is high-T c superconductivity in cuprate compounds © and the 
recently discovered iron-arsenic alloys ©. In order to capture the essential physics of such 
systems, the fermionic Hubbard Hamiltonian (3) has been introduced as a fundamental model 
describing interacting electrons in a periodic potential (UJ|4]). In a real solid, however, the effects 
of interest are typically complicated by e.g. multiple bands and orbital degrees of freedom, im- 
purities and the long-range nature of Coulomb interactions which becomes especially relevant 
close to a metal to insulator transition. It is therefore crucial to probe this fundamental model 
Hamiltonian in a controllable and clean experimental setting. Ultracold atoms in optical lattices 
provide such a defect-free system (|5]|6]), in which the relevant parameters can be independently 
controlled, allowing quantitative comparisons of the experiment with modern quantum many- 
body theories. For the case of bosonic particles a series of experiments carried out in the 
regime of the superfluid to Mott insulator transition (191 17011771) have demonstrated the versatility 
of ultracold quantum gases in this respect. For both bosonic and fermionic systems, the entrance 
into a Mott insulating state is signaled by a vanishing compressibility, which can in principle 
be probed experimentally by testing the response of the system to a change in external confine- 
ment. This is probably the most straightforward way to identify the interaction-induced Mott 
insulator and to distinguish it, e.g., from a disorder induced Anderson insulator (T721I731I741) . 
In a solid, however, the corresponding compressibility can usually not be measured directly, 
since a compression of the crystalline lattice by an external force does not change the number 
of electrons per unit cell (see also SOM.l). 

In this work, non-interacting and repulsively interacting spin mixtures of fermionic atoms 
deep in the degenerate regime are studied in a three-dimensional optical lattice. In the experi- 
ment, we are able to independently vary the interaction strength between the fermions using a 
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Feshbach resonance, as well as the lattice depth and the external harmonic confinement of the 
quantum gas. By monitoring the in-trap density distribution of the fermionic atoms for increas- 
ing harmonic confinements, we directly probe the compressibility of the many-body system. 
This measurement allows us to clearly distinguish compressible metallic phases from globally 
incompressible states and reveals the strong influence of interactions on the density distribu- 
tion. For non-interacting clouds the system changes continuously from a purely metallic state 
into a globally incompressible band-insulating state with increasing confinement. For repul- 
sive interactions, we find the cloud size to be significantly larger than in the non-interacting 
case, indicating the resistance of the system to compression. For strong repulsion, the system 
evolves from a metallic state into a Mott insulating state and eventually a band insulator as the 
compression increases. 

In order to probe the local on-site physics of the system, we measure the fraction of atoms on 
doubly occupied lattice sites for different experimental parameters. For medium compressions 
and strong repulsive interactions, we find an almost order of magnitude suppression of the frac- 
tion of doubly occupied sites compared to the non-interacting case. For strong compressions 
the systems enters a band-insulating state and the fraction of atoms on doubly occupied sites be- 
comes comparable for all interactions. In previous experiments, a suppression of the number of 
doubly occupied sites was demonstrated for increasing interaction strength for bosons (1751) and 
fermions (16) at fixed harmonic confinement, signaling the entrance into a strongly interacting 
regime. 

The experimentally observed density distributions and fractions of doubly occupied sites 
are compared to numerical calculations using Dynamical Mean Field Theory (DMFT) (|77l 
I7ffll79ll20|) . DMFT is a central method of solid state theory and is widely used to obtain ab- 
initio descriptions of strongly correlated materials (|7#1) . This comparison of DMFT predictions 
with experiments on ultracold fermions in optical lattices constitutes the first parameter-free 
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experimental test of the validity of DMFT in a three-dimensional system. 

Theoretical model. Restricting our discussion to the lowest energy band of a simple cu- 
bic 3D optical lattice, the fermionic quantum gas mixture can be modeled via the Hubbard- 
Hamiltonian (3) together with an additional term describing the potential energy due to the 
underlying harmonic potential: 

H = -J + u *ki"i,T 

+ V t J2(% + i l + r r 2 il)(nid + ^)- 

i 

Here the indices i,j denote different lattice sites in the three-dimensional system (i = 
iy, iz)), (i, j) neighboring lattice sites, a 6 {|, T} me two different spin states, J the tunnel- 
ing matrix element and U the effective on-site interaction. The operators Ci ta (c\ a ) correspond 
to the annihilation (creation) operators of a fermion in spin state a on the iih. lattice site and n i;(T 
counts the number of corresponding atoms on the ith lattice site. The strength of the harmonic 
confinement is parameterized by the energy offset between two adjacent lattice sites at the trap 
center V t = \muj\d 2 , with co± = uj x = cu y uj z being the horizontal trap frequency and d the 
lattice constant. The constant aspect ratio of the trap is denoted by 7 = lu z /lu±. Due to the Pauli 
principle every lattice site can be occupied by at most one atom per spin state. 

The quantum phases of the Hubbard model with harmonic confinement are governed by 
the interplay between three energy scales: kinetic energy, whose scale is given by the lat- 
tice bandwidth 12 J in three dimensions, interaction energy U, and the strength of the har- 
monic confinement, which can conveniently be expressed by the characteristic trap energy 
Et = VJ(7iV (T /(47r/3)) 2 / 3 , which denotes the Fermi-energy of a non-interacting cloud in the 
zero-tunneling limit with N a being the number of atoms per spin state (JVj, = N^). The char- 
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acteristic trap energy depends both on atom number and trap frequency via E t oc cu 2 , Na' and 
describes the effective compression of the quantum gas, which is controlled in the experiment 
by changing the trapping potential. 

Depending on which term in the Hamiltonian dominates, different kinds of many-body 
ground states can occur in the trap center (see Fig.Q]). For the case of weak interactions in a shal- 
low trap U <^E t <^ 12 J the Fermi energy is smaller than the lattice bandwidth {Ep < 12 J) and 
the atoms are delocalized in order to minimize their kinetic energy. This leads to compressible 
metallic states with central filling n 0;a < 1 (FigHK), where the local filling factor n it<T = (h i;(T ) 
denotes the average occupation per spin state of a given lattice site. A dominating repulsive in- 
teraction U ^>12 J and U^$> E t suppresses the double occupation of lattice sites and can lead to 
either Fermi-liquid (n 0i(7 < 1/2) or Mott-insulating (n 0:O - = 1/2) states in the center of the trap 
(FigOjB), depending on the ratio of kinetic to characteristic trap energy. Stronger compressions 
lead to higher filling factors, ultimately (E t S>12 J, E t ^$>U) resulting in an incompressible band 
insulator with unity central filling at T = (Fig[TJC). 

Finite temperature reduces all filling factors and enlarges the cloud size, as the system needs 
to accommodate the corresponding entropy. Furthermore, in the trap the filling always varies 
smoothly from a maximum at the trap center to zero at the edges of the cloud. For a dominating 
trap and strong repulsive interaction at low temperature (E t > U > 12 J), the interplay between 
the different terms in the Hamiltonian gives rise to a wedding-cake like structure (see also 
Fig. |4ji, F) consisting of a band-insulating core (n 0i(T ~ 1) surrounded by a metallic shell (1/2 < 
rtj iCT < 1), a Mott-insulating shell (n^ = 1/2) and a further metallic shell (rij jCT < 1/2) (19). The 
outermost shell remains always metallic, independent of interaction and confinement, only its 
thickness varies (see also SOM.8). 
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Experimental setup. Our experiments use an equal mixture of quantum degenerate fermionic 
40 K atoms in the two hyperfine states \F } m F ) = ||, — |) = ||) and ||, — |) = |f). By over- 
lapping two horizontal, elliptical laser beams (A = 1030 nm) a pancake-shaped dipole trap with 
an aspect ratio 7 w 4 is formed. Using evaporative cooling in this trap, we reach temperatures 
of T/Tp = 0.15(3) with 1.5 — 2.5 x 10 5 potassium atoms. A Feshbach resonance located 
at 202.1 G (|2~7T) is used to tune the scattering length between the two spin states and thereby 
control the on-site interaction U . The creation of the spin mixture and the last evaporation step 
are performed either above the resonance (220 G), giving access to non-interacting (209.9 G) as 
well as repulsively interacting clouds with a < 150 a (B < 260 G), or below the resonance 
(165 G), where larger scattering lengths up to a = 300 a (191.3 G) can be reached. We have 
found that a further approach to the Feshbach resonance is hindered by enhanced losses and 
heating in the lattice (see (|22l ). 

After evaporation, the dipole trap depth is again slightly increased and the magnetic field is 
tuned to the desired value. Subsequently, a blue detuned 3D optical lattice (A; at = 738 nm) with 
simple cubic symmetry is linearly increased within 10 ms to a potential depth of V[ at = 1 E r . 
Here E r = h 2 /{2m\ 2 at ) denotes the recoil energy, which sets a natural energy scale for the 
depth of the optical lattice © and m is the mass of a single atom. The combination of a red- 
detuned optical dipole trap and a blue-detuned lattice potential allows us to vary lattice depth 
and external confinement independently. Thereby, the compression of the atomic cloud can be 
adjusted over a wide parameter range and horizontal trap frequencies as low as lu± ~ 2n x 20 Hz 
can be reached in the presence of the lattice, thus allowing for metallic states with high atom 
numbers. To monitor the in- situ density distribution under different external confinements, 
we then ramp the dipole trap depth in 100 ms to the desired external harmonic confinement 
(tu± = 2n x 20 — 120 Hz), followed by a linear increase of the optical lattice depth to Vi at = 8 E r 
in 50 ms. 
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An in-situ image of the cloud is subsequently taken along the short axis of the trap using 
phase-contrast imaging (l23l) at detunings of A = 2ix x (200 — 330) MHz after a hold time of 
12 ms in the lattice. From this picture the cloud size R = (rj_) is extracted using adapted 2D 
Fermi fits (see Supporting Online Material (SOM.5)). 

As phase-contrast imaging is non-destructive, it allows us to also measure the quasi-mo- 
mentum distribution of the atoms in the same experimental run using a band-mapping tech- 
nique ([2"?ll2"5ll2"oT) . To this end, the lattice is ramped down in 200 /is and a standard absorption 
image is taken after 10 ms time- of- flight. 

All experimental data are compared to numerical calculations, in which the DMFT equations 
of the homogeneous model are solved for a wide range of temperatures and chemical potentials 
using a numerical renormalization group approach (|27][2#|) (see SOM.7 for details). As shown 
in (|29ll79l) . the trapped system can be approximated to very high accuracy by the uniform 
system through a local density approximation (LDA) even close to the boundary between metal 
and insulator. For a comparison with the experimental results it is convenient to express the 
cloud size R in rescaled units R sc = i?/(7iV cr ) 1 / 3 , along with the dimensionless compression 
E t /\2J. In these units, the cloud size depends only on the interaction strength U/12J and the 
entropy. In all calculations, we use the entropy determined from a non-interacting Fermi gas in 
an harmonic trap at an initial temperature T/Tp and assume adiabatic lattice loading 

A comparison between the numerically calculated density distributions, the corresponding 
column densities and the experimentally measured ones is presented in Figure [2] for differ- 
ent confinements and interactions at constant atom number. At low compression (first row) 
all shown distributions are purely metallic, as the central filling of n 0)Cr = 0.65 in the non- 
interacting case (black curve) is reduced by repulsive interactions to below half filling (blue, 
red). Due to the fixed atom number, this reduction in density causes an enlargement of the 
cloud, which is clearly visible in the experiment. In the second row the compression is slightly 
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higher (E t /YlJ = 0.35) and a Mott-insulating core with half filling starts to form for repulsive 
interactions (A2), while the non-interacting distribution remains metallic. For high compres- 
sions (E t /12J = 1.4, last row) a band-insulating core has formed in the non-interacting case 
(A3), while the repulsively interacting clouds show a metallic core. For strong interactions 
(red) the remains of a Mott-insulating ring can be seen (A3). All experimental profiles agree 
well with theory, small deviations at low compressions could be attributed to an anharmonicity 
of the confining potential for very shallow traps. 

This general behavior can be quantified further by plotting the cloud size R sc (dots) in 
rescaled units as a function of the characteristic trap energy E t oc cu \Nl /z (see Fig.©. Addi- 
tionally, the global compressibility k Rsc = — ^- df ^f/i2j) ( see Fig.|4]and SOM.l) of the system 
can be extracted from these measurements using linear fits to four consecutive data points to 
determine the derivative. In the non-interacting case we find the cloud sizes (Fig. [3l black dots) 
to decrease continuously with compression until the characteristic trap energy roughly equals 
the lattice bandwidth (E t /12J ~ 1). For stronger confinements the compressibility approaches 
zero (Fig. UK), as almost all atoms are in the band insulating regime while the surrounding 
metallic shell becomes negligible. The corresponding quasi-momentum distribution (Fig. [3}\- 
E) changes gradually from a partially filled first Brillouin zone, characteristic for a metal, to 
an almost evenly filled first Brillouin zone for increasing compressions, as expected for a band 
insulator. We would like to note, however, that a band-mapping technique reveals only the 
relative occupations of the extended Bloch states. For an inhomogeneous system it therefore 
yields no information about the real-space density and especially cannot distinguish insulating 
from compressible states, e.g. non-equilibrium states in which the atomic wavefunctions are 
localized to single lattice sites. In contrast, the measurements shown here directly demonstrate 
the global incompressibility of the fermionic band insulator, in excellent agreement with the 
theoretical expectation for a non-interacting Fermi gas (black line). 
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The green, blue and red dots in figures [3] and @] represent the size and compressibility of 
repulsively interacting clouds with £7/12 J = 0.5,1 and 1.5 in comparison with the DMFT 
calculations (lines). For moderately repulsive interaction (U/12J = 0.5,1) the cloud size 
is clearly bigger than in the non-interacting case but eventually reaches the size of the band 
insulator. For stronger repulsive interactions (U/12J = 1.5) we find the onset of a region 
(0.5 < E t jYlJ < 0.7) where the cloud size decreases only slightly with increasing harmonic 
confinement, denoting a very small compressibility, whereas for stronger confinements the 
compressibility increases again. This is consistent with the formation of an incompressible 
Mott-insulating core with half filling in the center of the trap, surrounded by a compressible 
metallic shell, as can be seen in the corresponding in-trap density profiles (see Fig. HE, F). For 
higher confinements an additional metallic core (1/2 < n i a < 1) starts to form in the center of 
the trap. A local minimum in the global compressibility is in fact a genuine characteristic of 
a Mott-insulator and for large U and low temperature, we expect the global compressibility in 
the middle of the Mott region to vanish as 1/U 2 (see SOM.l). The experimental data, indeed, 
show an indication of this behavior (see Fig.|4]C) for increasing interactions. For i? 4 /12J~0.5 
a minimum in the compressibility is observed, followed by an increase of the compressibility 
around _E t /12J~0.8, slightly earlier than predicted by theory. 

When the system is compressed even further all cloud sizes approach that of a band insu- 
lating state and all compressibilities tend to zero. As can be seen in the theory predictions, 
for strong confinement the repulsively interacting clouds can, however, become slightly smaller 
than the non-interacting one due to Pomeranchuk cooling (30). At the same average entropy per 
particle, the interacting system has a considerably lower temperature in the lattice, as the spin 
entropy is enhanced due to interactions (see SOM.3). In the experiment this feature is barely 
visible, as a second effect becomes important: At very high compressions (E t /12J > 2) the 
second Bloch band gets slightly populated during the lattice ramp up, which leads to smaller 
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cloud sizes for all interactions, as a small number of atoms in a nearly empty band can carry a 
considerable amount of entropy. 

Overall, we find the measured cloud sizes to be in very good quantitative agreement with the 
theoretical calculations up to U/12J = 1.5 (B = 175G). Nevertheless, for repulsive interactions 
and medium compression (E t /12J « 0.5) the cloud size is slightly bigger than the theoretical 
expectation. The discrepancies become more prominent for stronger interactions, i.e. when 
tuning the scattering length to even more positive values below the Feshbach resonance. This 
could be caused by non-equilibrium dynamics in the formation of a Mott-insulating state for 
strong interactions (TJ/l) or may be an effect not covered by the simple single band Hubbard 
model or the DMFT calculations and requires further investigation. 

In order to ensure that the used lattice loading time of 50 ms is adiabatic, we have measured 
the resulting in-situ cloud size as a function of ramp time (see inset Fig. I3F) in the regime 
around E t /l2J 0.5, where the differences between experiment and theory are biggest. In 
this regime the cloud size shrinks during the loading of the lattice, a too fast loading thereby 
results in a larger cloud. We find the cloud size to decrease significantly for ramp times up 
to 20 ms after which no further decrease is observed, which indicates adiabaticity for the used 
ramp time. However, a second longer timescale, which could become more relevant for stronger 
interactions (1321 . cannot be ruled out. In addition, the temperatures before loading into the 
lattice and after a return to the dipole trap with a reversed sequence are compared. We find 
a rise in temperature between 0.010(5) T/Tp for a non-interacting cloud and 0.05(2) T/Tp 
for a medium repulsion of U/12J = 1 at compressions around E t /12J w 0.5. Even for 
the highest compressions the temperature increase stays below 0.11(2) T/Tp. Taking these 
temperature increases into account, the good agreement between the experimental data and the 
numerical calculations, which assume adiabatic loading and an initial temperature of T/Tp = 
0.15, indicates that our actual initial temperatures lie rather at the lower end of the measured 



10 



temperature range T/T F = 0.15(3). 

The theoretical calculations for the compressibility shown in Fig. 0J3 demonstrate that the 
entrance into a Mott insulating state, signaled by a minimum in the global compressibility, 
happens rather rapidly for initial temperatures in the range of 0.15 < T/T F < 0.2. At these 
temperatures, the entropy per particle is larger than k B 2 In 2, which is even larger than the max- 
imum possible entropy per particle of a half-filled Hubbard model in the homogeneous system 
(reached only at temperatures in the lattice ksTiat 3> U). In the trap, however, a large frac- 
tion of the entropy is accumulated in the metallic shells at the edges of the atomic cloud where 
the diluted atoms have a large configurational entropy (see SOM.4). Therefore the temperature 
remains on the order of ksT^t ~ J <C U in the Mott insulating regime (see SOM.3). 

In addition to the global compressibility measurements, the fraction of atoms on doubly 
occupied lattice sites (pair fraction p) has been measured for magnetic fields above the Feshbach 
resonance (U/12J = 0, 0.5, 1) by converting all atoms on doubly occupied sites into molecules 
using a magnetic field ramp (0.2ms/G) over the Feshbach resonance (|27][76]). In this case, the 
lattice depth is further increased to 20 E r in 200 /is to prevent tunneling of the atoms during 
the field ramp. After the ramp, the lattice depth is linearly decreased to zero in 200 fis and 
the number of remaining atoms is recorded by time-of-flight absorption imaging, yielding the 
number of singly occupied lattice sites. The difference in atom number with and without the 
magnetic field ramp normalized to the total atom number gives the desired fraction p, which is 
plotted in Fig. [5K including corrections which account for atom losses during the measurement 
sequence (see also SOM.6). The fraction of atoms on doubly occupied sites gives insight into 
the local on-site physics of the system. In combination with the in-situ size measurements, this 
fraction can be compared for different interaction strengths at constant average density of the 
system. The average density n a = (7 N a /^irR 3 ) x l/n BI , where n BI denotes the average 
density of a pure band-insulator at T = 0, can be calculated from the recorded cloud size R and 
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is shown in Fig. 03. 

In the limit of very low compression (weak confinement) the average density is small 
n a < 0.1 and p tends to zero. For intermediate densities the fraction of doubly occupied 
sites depends crucially on the interaction. At a constant mean density of n a = 0.4, we find a 
pair fraction of 40% for a non-interacting cloud and around 5% for a medium repulsive inter- 
action U/12J = 1 (see yellow circles in Fig.0\, B). In this regime of repulsive interactions at 
medium compression it is energetically favorable to reduce the number of doubly occupied sites 
although this costs potential and kinetic energy. As a consequence, different compressions are 
needed to reach the same average density for different interactions, as can be seen in Fig. 03. 

For strong compressions, the measured pair fraction becomes comparable for all interactions 
(see Fig. 0, as the potential energy in the trap becomes larger compared to U and all atom 
distributions are expected to contain a large band-insulating core (see SOM.8). Ultimately, the 
pair fraction (average density) is limited to values smaller than 60% (n a < 0.7) due to the finite 
entropy per particle, which reduces the filling factor in the band-insulating state. While the non- 
interacting and slightly repulsively interacting curve U/ 12 J — 0.5 match the DMFT results for 
an inital temperature of T/Tp = 0.15, we see deviations for stronger repulsive interactions 
(U/12J = 1). In this case the measured pair fraction is in general ~ 10% higher than predicted 
by theory, nevertheless the qualitative behavior agrees very well. We note that a suppressed pair 
fraction in comparison with the non-interacting case occurs for all temperatures in the lattice 
below &b Tiat ~ U, regardless of the formation of an incompressible Mott insulating phase in 
the inhomogeneous system (see SOM.2). Furthermore, the pair fraction vanishes even for a 
compressible purely metallic phase with n iy(T < 1/2 in the strongly interacting regime. 

Summary and Outlook. In conclusion, we have investigated a spin mixture of repulsively 
interacting fermionic atoms in a 3D optical lattice. Using a novel measurement technique, we 
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have been able to directly determine the global compressibility of the many-body quantum sys- 
tem and have explored the different regimes of the interacting mixture from a Fermi liquid to a 
Mott and band insulating state upon increasing harmonic confinements and increasing interac- 
tions. Additionally, the fraction of atoms on doubly occupied sites has been recorded in order to 
determine the local on-site behavior of the system for a large range of compressions. Our mea- 
surements are in good agreement with the results predicted by DMFT for interaction strengths 
up to U /12J = 1.5 and initial temperatures of the mixture of T/Tp = 0.15. For stronger in- 
teractions, we find deviations from the expected cloud radius in the low compression regime, 
which could point to interesting non-equilibrium dynamics in the loading process (TJTT) . that 
hinders approaching the expected equilibrium many-body state of the system in this parameter 
regime. Our numerical calculations also show that low initial temperatures of T/T F < 0.15 are 
needed to enter an incompressible many-body Mott insulating state, whereas a suppressed pair 
fraction compared to the non-interacting state can already be observed at much higher temper- 
atures and for metallic states. This behavior is in fact similar to the results obtained in recent 
calculations and experiments on the melting of incompressible bosonic Mott insulating shells 
for increasing temperatures (13JIIJ4[) . 

Our measurements are an important step in the direction of analyzing fermionic many-body 
systems with repulsive interactions in a lattice. For initial entropies lower than S/N < ks In 2, 
one expects the system to enter an antiferromagnetically ordered phase, as the temperature of 
the quantum gas can then drop below the superexchange coupling that mediates an effective 
spin-spin interaction between the particles (T2Q][30][35ll3oT) . This would open the path to the 
investigation of quantum magnetism with ultracold atoms (TJTI) . being an encouraging starting 
point to ultimately determine the low-temperature phase diagram of the Hubbard model (17113^). 
This includes the search for a c?-wave superconducting phase (|39|) that is believed to emerge 
from within the two-dimensional Hubbard model. 
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Figure 1: Relevant phases of the Hubbard model with an inhomogeneous trapping potential 
for a spin mixture at T = 0. A schematic representation is shown in the left column. The 
center column displays the corresponding in-trap density profiles and the right column outlines 
the distribution of singly and doubly occupied lattice sites after a rapid projection into the zero 
tunneling limit J = 0, with p denoting the total fraction of atoms on doubly occupied lattice 
sites. 
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Figure 2: Comparison of in-trap density profiles. Calculated radial density profiles for different 
compressions (harmonic confinements) E t /\2J (left column (A1-A3)), corresponding column 
densities obtained after integration over the z— axis and convolution with the point spread func- 
tion of our imaging system (center column (B1-B3)) and experimental results (azimuthally av- 
eraged over > 5 shots) for three different interaction strengths U/12J = (black), U /12J = 1 
(blue) and U/12J = 1.5 (red) (C1-C3). At small compressions (A, B) the calculated den- 
sity profiles for U /12J = 1 and U /12J = 1.5 are indistinguishable, as in both cases double 
occupations are almost completely suppressed for n, >cr < 1/2. 
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Figure 3: Cloud sizes of the interacting spin mixture versus compression. Measured cloud size 
R sc in a Vi at = 8 E r deep lattice as a function of the external trapping potential for various 
interactions U/12J = (black), U/12J = 0.5 (green), U/12J = 1 (blue), U/12J = 1.5 
(red)). Dots denote single experimental shots, lines the theoretical expectation from DMFT for 
T/Tp = 0.15 prior to loading. The insets (A-E) show the quasi-momentum distribution of 
the non-interacting clouds (averaged over several shots). (F) Resulting cloud size for different 
lattice ramp times at E t /12J = 0.4 for a non-interacting and an interacting Fermi gas. The 
arrow marks the ramp time of 50 ms used in the experiment. 
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Figure 4: Compressibility and in-trap density distribution. (A-C) Global compressibility k Rsc 
of the atom cloud for various interactions. Dots denote the result of linear fits on the measured 
data with the error bars being the fit uncertainty for three interaction strengths ((A) U/12 J — 0, 
(B) U /12J = 1, (C) U /12J = 1.5) Solid lines display the theoretically expected results for an 
initial temperature of T/T F = 0.15. The influence of the initial temperature on the calculated 
compressibility is shown in D for U/12 J = 1.5. The corresponding density distributions are 
plotted in E, F with r denoting the distance to the trap center (see also SOM.8). The red lines 
mark the region where a Mott-insulating core has formed in the center of the trap and the global 
compressibility is reduced. 



17 



60 




Compression E t H2J 

Figure 5: Fraction of atoms on doubly occupied sites versus compression for different interac- 
tion strengths (black U/12J = 0, green U/12J = 0.5, blue U/12J = 1) (A). The inset (B) 
displays the evolution of the average density for increasing compression and different interac- 
tion strengths (as in A, red U/12J = 1.5), allowing one to compare the fraction of doubly 
occupied sites at constant average density. The yellow circles denoted by a and a' indicate 
the fraction of atoms in doubly occupied sites for a constant average density of n a = 0.4 and 
U/12J = and U/12J = 1. 
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Metallic and Insulating Phases of Repulsively 
Interacting Fermions in a 3D Optical Lattice 

Supporting online materials 

1 Compressibility of atom clouds 

The concept of compressibility is used in rather different ways in physics. Often the term 
refers to the change of volume when pressure is applied, k v = — To investigate the 

effects of strong interactions on various phases, a different concept, usually called electronic 
compressibility, is very useful. Here one tracks the change of the number of electrons per unit 
cell n of a solid as a function of the chemical potential, n n = ^f^. A central property of a 
Mott insulator is its incompressibility: n n vanishes exponentially for low temperatures. The 
compressibility allows one to distinguish a Mott insulator, e.g., from an Anderson insulator of 
non-interacting particles in the presence of strong disorder. n n can be measured only indirectly 
in three-dimensional solids, e.g. by comparing shifts in photoemission spectra (SI) for samples 
with different doping. The interpretation of such an experiment is, however, difficult, as a 
doping of a sample also changes other contributions to the total energy (e.g. the Coulomb 
energy) which also lead to a shift of the spectra. A more direct measurement of K n can be 
obtained in thin films or field effect transistors where the density of electrons can be changed 
using a back gate, see e.g. (S2). 

It is important to realize that in solids the compressibility of the volume k v is usually unre- 
lated to the electronic compressibility n n as the number of electrons per unit cell is not modified 
when pressure is applied. The situation is completely different for atomic clouds. When an 
atomic cloud is compressed by a change of the trapping potential V t the density of atoms does 
change considerably as the radius of the cloud R will be reduced. Therefore the compression 
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of an atom cloud gives valuable information on the properties of the various phases realized in 
the trap. 

The global compressibility of the cloud can be defined as = —-gg^jr, where R = 
J r 2 n c (r)d 2 r /N is the typical radius of the cloud (N = 2N a ) and n c (r) is the column den- 
sity. For simplicity, we use an isotropic trap (7 = 1) in the following. The powers of R in 
the prefactors are chosen such that the quantities are well-defined in the thermodynamic limit 
N a — > 00 for fixed E t . In rescaled units the corresponding definition is k Rsc = —^- Q^f/nj) = 
12J(|tt) 2 / 3 ^ 

As has been shown in ( S3,S4), the local density of the system is well described by the density 
of a homogeneous system with the chemical potential /i — V t r 2 (LDA approximation). When V t 
is changed, both the chemical potential n and the temperature T are modified. Using that both 
particle number and entropy remain constant, one can eliminate d/j/dV t and dT/ dV t to express 
the total compressibility of the trap in terms of grand-canonical derivatives of the density n and 
entropy s per lattice site 

where r is a typical radius defined by 

(rl - r 2 ? ( ^ 2 

- 2 + 

r 



< = < + —7 (2) 

2 / 8N 8S _ 8N 8S 



8N 

c \dfj, dT dT d/j, 

with r 2 c = Jr 2 ^d 3 r/ J ^d 3 r and r 2 = fr 2 j^d 3 r/f ^d 3 r and the total entropy S = 
J sd 3 r. For low temperature one can neglect the rise of the temperature when V t is increased 
and one obtains r ~ r c . Under the experimental conditions, however, one has to take the 
heating effects into account to describe the experiments correctly. 

Eq. [T] can be used to estimate the size of k,r deep in the Mott insulating phase when the 
Mott insulating core is surrounded by a compressible shell of width Ar. The metallic shell is 
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approximately located at r and its width is inversely proportional to the slope of the potential, 
Ar ~ J / (VoTo). Due to particle number conservation reflected in the r 2 , counterterm in Eq. (0Q), 
k r is proportional to (Ar) 2 , k r ~ ^ Ar J r ° . Using that V rQ ~ [/ and r ~ i? ~ iV 1 / 3 we 
obtain the expected result that n R vanishes rapidly for large U, k r ~ reaching J/U 2 in 

the center of the Mott plateau, see Figure SI. When V t is further increased, particles from the 
outer shell can move to the center of the trap and the counter-term r 2 becomes much smaller. 
In this regime n R is linear in Ar and therefore n R ~ l/U until the center of the trap is filled 
with doubly occupied sites. In the low density limit, V t — ► 0, k r diverges proportional to 1/ \fV t 
at low temperatures as interactions can be neglected and R ~ N 1 l' i (E t /12J)~ l / A . The above 
described behavior of k r is shown in Figure SI. 



(J Ef" 




12 J L//2 (J 



5 

Figure SI: kj? as a function of i?t. The qualitative behavior is shown using the DMFT results for 
U/12J = 2.5 for an initial temperature of T/Tp = 0.15 as an example. The distinct minimum 
in the global compressibility k r is a clear signature of the incompressible Mott-insulator. For 
lower temperatures, the minima and maxima become more pronounced, see Fig.4 of the main 
paper. 
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E t ii2j 

Figure S2: Calculated global compressibility k Rsc in rescaled units as a function of E t /12J for 
different interaction strengths (U/12J = . . . 2.5). All curves are calculated for the same initial 
temperature T/T F = 0.15 prior to loading of the lattice. 

Figure S2 shows the calculated global compressibility for various interactions and an initial 
temperature of T/Tp = 0.15. Starting at U/12J = 1 a minimum in the global compressibil- 
ity appears, as an incompressible Mott-insulator forms in the center of trap (for the uniform 
system at T = critical interaction for the metal to insulator transition, within DMFT, is 
U/12J = 1.26 (S4J). For larger interactions the minimum becomes broader and deeper as 
the Mott-insulting core increases in size (see density profiles in SOM.8). The remaining finite 
global compressibility is dominated by the compressible metallic shell around the incompress- 
ible Mott-insulator and decreases like oc 1/U 2 . 
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2 Melting of the Mott-insulator for increasing temperature 




0.2 0.4 0.6 0.8 1.0 

E t h2J 



Figure S3: Melting of the Mott insulator with U/12J = 1.5 for increasing temperature. A: 
calculated central density, B: corresponding pair fraction p. 

In Figure S3A the central in-trap density is plotted for various initial temperatures as a function 
of compression. The vanishing slope at E t /12J rs 0.5 for low enough temperatures directly 
shows the local incompressibility (k„ rs 0) of the Mott insulator, as the chemical potential 
increases with higher compressions. For temperatures around T/Tp « 0.15 the Mott insulator 
melts and the local compressibility strongly increases. The pair fraction on the other hand, 
which is shown in Figure S3B remains strongly suppressed and changes by less than 1% over 
the whole temperature range, showing no signature of the crossover into the Mott insulator. 
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3 Pomeranchuk cooling 




Figure S4: Calculated temperature and cloud size in the lattice. A: Calculated temperature 
Ti at / 12 J in the lattice for various interaction strengths at a constant initial temperature of 
T/Tp = 0.15 prior to the loading. Small wiggles in the curves are artifacts of the interpola- 
tion of the entropy, see SOM.7. B: Calculated cloud size as a function of compression E t /Y2J 
for the same interactions as in A. At compressions larger than one the size of the repulsively 
interacting clouds becomes smaller than in the noninteracting cloud due to the Pomeranchuk 
effect. 



The final temperature in the lattice, which is plotted in Figure S4A in rescaled units, is a function 
of the initial entropy per particle, the compression and the interaction. Due to the Pomeranchuk 
effect, the temperature in the lattice T[ at decreases for stronger interactions as the particles 
become more localized, which increases the spin entropy. Therefore the cloud size at high 
compressions can become slightly smaller for medium interactions than for zero interaction, as 
can be seen in Figure S4B. 
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4 Entropy distribution 




Distance from trap center r (d) 

Figure S5: Calculated entropy distributions in the lattice. Solid (dashed) lines show the entropy 
per lattice site (per particle) for initial temperatures of T/Tp = 0.07 (black) and T/Tp = 0.15 
(red) and strong repulsive interaction U/12J = 1.5. 

In Figure S5A the entropy distribution of a purely metallic state with less than half filling in the 
center of the trap is shown. At higher compression (B) a Mott-insulator with unity filling and 
k B In 2 entropy per particle has formed in the center of the trap even in the case of T / T F = 0.15, 
for which the average entropy per particle is above fc B 2 hi 2. This is possible only due to the 
inhomogeneity of the system, as most of the entropy is carried by the metallic shells where the 
entropy per particle can diverge. For high compressions (D) a band-insulating core has formed 
and for low enough temperatures (black) nearly all entropy is carried by the surrounding shells. 
The small dip at r = 20 is a remnant of the Mott-insulating shell between the two metallic 
shells. Note that entropy per lattice site and entropy per particle are equal at half filling. 
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5 Fitting procedure 



To extract the cloud size R from the phase contrast images, we fit the 2D image of the cloud 
with the following adapted Fermi-Fit function: 

-100e ^ +b + cJ [ J c) + [V J c) 

with Li 2 being the di-logarifhm and x c , y c , a x , a y , a, b, c free fit parameters. The last term mod- 
els a broad funnel-shaped background, which is an artifact of phase-contrast imaging. As can 
be seen in Figure S6, this function describes our measured curves much better than a Gaus- 
sian. Indeed the adapted Fermi-Fits yield on average 8% (non-interacting cloud) and 23% 
(interacting cloud with U/12J = 1.5) smaller squared residuals compared to a Gaussian fit 
function including the last background term. Using this fit-function, the imaged cloud size 
R = a/ (rj_) after integration over the propagation axis of the imaging laser beam is given by 
R = y^l.264 2 (cr^ + cr^) — rj 2 , where rj denotes the width of the point spread function of the 
imaging setup, which is well below one third of the smallest used cloud size. 
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Figure S6: Comparison of different fit functions. The azimuthally averaged data (blue dots) is 
shown together with a fitted Gaussian (red line) and an adapted Fermi-fit (black line). Both fits 
were performed on the full 2D data before averaging. The influence of the small deviation at 
the trap center on the resulting atom number (cloud size) is below 1% (0.1%), as the statistical 
weight of each averaged data point increases linearly with the distance r from the cloud center. 

6 Atom number calibration and fraction of atoms on doubly 
occupied sites 

The atom numbers in the manuscript were calibrated such that the theoretically calculated renor- 
malized cloud size for non-interacting atoms at high compressions was reproduced. This atom 
number calibration agrees within 10% with both the theoretical expectation when taking the 
polarization of the probe beam and the corresponding Clebsch-Gordan coefficients of the tran- 
sition used for phase contrast imaging into account and a second independent calibration, in 
which the cloud size of a non-interacting Fermi gas in a pure dipole trap was measured as a 
function of the trap frequency and compared to the theoretically expected scaling. 

When measuring the fraction of atoms on doubly occupied lattice sites, we need to take 
a strong reduction of the pair lifetime for high lattice depths Vi at = 20 — 25E r into account. 
During the 15 ms hold time used for the Feshbach ramps, doubly occupied sites can be lost 
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due to a light assisted collision process caused by the blue detuned lattice light (S5,S6). This 
process selectively affects only doubly occupied sites and therefore yields an underestimation 
of the measured pair fraction p, for which the experimental data has been accordingly corrected. 
The observed atom losses range from 15% for non-interacting atoms to below 5% for strong 
repulsive interaction and are consistent with the independently measured lifetimes of 70(10) ms 
(a = 0a ) resp. 310(70) ms (a = 150a ) for doubly occupied sites in a deep lattice V tat = 2AE r . 

7 Numerical implementation 

To obtain the theoretical results shown in the paper, we use the numerical renormalization 
group method (NRG) to solve the dynamical mean field equations for the Hubbard model in 
the thermodynamic limit, see Refs. (S7,S8). The NRG calculations were carried out using a 
logarithmic discretization of the conduction band with discretization parameter A/6 J = 1.5, 
N sh = 78 energy shells and 660 retained states per energy shell. From the solution we obtain 
the average number of atoms per site and spin, n(/x, T), and the average number of doubly 
occupied sites, n d (fj,, T), for a dense linear grid of chemical potentials and a logarithmic grid 
of temperatures. We interpolate the result. As the NRG code does not work very well for 
extremely low densities of particles or holes, we use the non-interacting result in this regime 
for n(n,T) (typically n < 0.03 for moderate T). It is, however, not possible to approximate 
n d (fi,T) by the non-interacting result in this regime. For finite temperature, the number of 
doubly occupied sites is proportional to n 2 for n — > with a proportionality factor which is in- 
dependent of T in this classical limit. We determine the proportionality factor c(U) numerically 
(c(U) « 0.39, 0.22, 0.13, 0.087, 0.064 for values of U ranging from U/12J = 0.5 to 2.5) and 
use n d w c(U)n 2 in the low-density limit. 

A direct calculation of the entropy within DMFT+NRG is subtle, therefore we use the 
thermodynamic relation dn/dT = ds/dfx to obtain s(fi,T) = J^dn/dT for /i < U/2 
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and s(/i, T) = s(U — fj,, T) for // > [7/2. The interpolation and integration errors due 
to this procedure are typically less than 5% as one can check by comparing s(fx = U/2) 
with In 2, the value expected in the Mott insulating phase. We then use the local density 
approximation (LDA), justified in (S3,S4)), to obtain the entropy per particle S/N(/j,,T) = 
J s(/i — V t r 2 , T)d 3 r/ J n(/2 — V t r 2 , T)d 3 r as a function of \i and T. Note that S/N does not 
depend on V t within LDA. For fixed V t , \x and S/N we then determine the temperature T and 
the particle number N = J n{ji — V t r 2 , T)d 3 r. Up to a trivial rescaling, all physical properties 
depend (within LDA) only on E t /(12J) oc V t N 2 ^. We can therefore directly obtain quantities 
like the fraction of atoms on doubly occupied sites p = N d /N = J n d (fi — V t r 2 , T)d 3 r/N, the 
rescaled column density n c (r)/A^ 1//3 = J n(fx — V t (r 2 + z 2 ),T)dz/N 1 ^ 3 or the effective radius 
of the cloud R/N 1 / 3 = yj J r 2 n c (r)d 2 r/ J n c (r)d 2 r as a function of E t /\2J from the above 
described procedure. 
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8 Density distributions 
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Figure S8: Calculated density distributions for different interactions and an initial temperature 
of T/Tp = 0.15. The small ridges for higher interactions are numerical artifacts. 



In the figures S7 and S8 the in-trap density distributions are plotted for various interaction 
strengths. The formation of the Mott-insulating core with half filling is clearly visible for both 
temperatures. 

References and Notes 

• SI. A. Ino, etal, Phys. Rev. Lett. 79, 2101 (1997). 

• S2. S. Hani, A. Yacoby, D. Mahalu, H. Shtrikman, Science 292, 1354 (2001). 

• S3. R. W. Helmes, T. A. Costi. A. Rosch, Phys. Rev. Lett. 100, 056403 (2008). 

• S4. R. W. Helmes, T. A. Costi, A. Rosch, Phys. Rev. Lett. 101, 066802 (2008). 

• S5. V. Vuletic, C. Chin, A. J. Kerman, S. Chu, Phys. Rev. Lett. 83, 943 (1999). 

• S6. P. S. Julienne, J. Res. Natl. Inst. Stand. Technol. 101, 487 (1996). 

• S7. R. Bulla, T. A. Costi, D. Vollhardt, Phys. Rev. B 64, 045103 (2001). 

• S8. R. Bulla, T. A. Costi, T. Pruschke, Rev. Mod. Phys. 80, 395 (2008). 



34 



